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Abstract 

We consider in details the effects of tlie 13-mixing (sin^^ia) and of the type of mass 
hierarchy/ordering (sign[Amf3]) on neutrino signals from the gravitational collapses 
of stars. The observables (characteristics of the energy spectra of fe and Vg, events) 
sensitive to sin^ ^13 and sign[Am^3] have been calculated. They include the ratio 
of average energies of the spectra, te = {E)/{E), the ratio of widths of the energy 
distributions, rr = T/f, the ratios of total numbers of i^e and events at low energies, 
S, and in the high energy tails, Rtaii- We construct and analyze scatter plots which 
show the predictions for the observables for different intervals of sin^ ^13 and signs 
of Ami2, taking into account uncertainties in the original neutrino spectra, the star 
density profile, etc.. Regions in the space of observables r^, rr, S, Rtaii exist in which 
certain mass hierarchy and intervals of sin^ ^13 can be identified or discriminated. 
We elaborate on the method of the high energy tails in the spectra of events. The 
conditions are formulated for which sin^ ^13 can be (i) measured, (ii) restricted from 
below, (iii) restricted from above. We comment on the possibility to determine sin^ ^13 
using the time dependence of the signals due to the propagation of the shock wave 
through the resonance layers of the star. We show that the appearance of the delayed 
Earth matter effect in one of the channels {ue or Ve) in combination with the undelayed 
effect in the other channel will allow to identify the shock wave appeareance and 
determine the mass hierarchy. 
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1 Introduction 



Collapsing stars (some of them appearing as supernova explosions) are the sources of neu- 
trinos of different flavors which can be used for oscillation/conversion experiments p. The 
structure of the neutrino mass spectrum and lepton mixing is imprinted into the detected 
signal (see |2]^C2] as an incomplete list of relevant works). Therefore, in principle, studying 
the properties of a supernova neutrino burst one can get information about 

• the values of parameters relevant for the solution of the solar neutrino problem, 

• the type of the mass hierarchy /ordering, 

• the 13-mixing parameter sin^6'i3, 

• the presence of sterile neutrinos, 

• new neutrino interactions. 

The first KamLAND results ^1] confirmed the LMA MSW as the dominant mechanism 
of the solar neutrino conversion. Further KamLAND measurements and solar neutrino 
studies will determine the corresponding oscillation parameters with rather good accuracy 
[TH t \W \ IT7j. This confirmation implies, in particular, that already in 1987, with the detec- 
tion of neutrinos from the supernova SN1987A, significant conversion effects were observed 
on supernova neutrinos JH]^|2ni- The identification of the neutrino mass hierarchy and 
the determination of 13-mixing (sin^ ^13) have become the main issues of further studies. 
Searches for sterile neutrinos and new neutrino interactions are also on the agenda. 

In this paper we will concentrate on the related subjects of the mass hierarchy and 
sin^6'i3. We will consider a three neutrinos system, assuming that sterile neutrinos, if they 
exist, produce negligible effect. 

The possibilities to study neutrino conversion effects using supernovae are wide, but not 
exempt of problems. The main difficulty originates from the fact that oscillation effects 
are proportional to the difference of the electron and non-electron neutrino fluxes originally 
produced inside the star, which are poorly known at the moment. The features of these 
fluxes depend on many details of the neutrino transport inside the star and, in general, on 
the type of progenitor star. 

There are two approaches to resolve the problem: 

1). perform a global fit of the data, determining both oscillation parameters and the param- 
eters of the original fluxes simultaneously. However, the number of unknown parameters 
which describe the energy spectra of the emitted neutrinos (temperatures, luminosities, 
pinching parameters) is rather large, and moreover, these quantities change with time dur- 
ing the burst. 
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Also degeneracies of parameters exist, so that variations of the oscillation parameters 
and of the parameters of the fluxes can produce the same observable effect. 

2). perform a (supernova) model independent analysis relying on some generic (and model 
independent) qualitative features of the fluxes. Namely: 

- the inequality of average energies (temperatures) of the original fluxes of neutrinos of 
different flavors; 

- the dominance of the electron neutrino flux in the initial phase of the burst (neutronization 
peak) 

- the pinching of the energy spectra 

- the approximate equality of the original u^, u^, u^-, fluxes. 

The first approach has been used recently in piU ITT]. while the study of ratios of numbers 
of events in specific energy intervals has been suggested in [HI CB] • 

In this paper we elaborate the second type of approach, suggesting specific methods and 
taking into account all the possible uncertainties. Both analytical and numerical analyses 
are performed. The paper is organized as follows. In sect. 2 we summarize our knowledge 
of the original neutrino spectra and of the matter distribution along the trajectory of the 
neutrinos in the star. We also introduce the scheme of neutrino masses and mixings. In 
sect. 3. the dependence of the conversion probabilities of supernova neutrinos on the mass 
hierarchy and the 13-mixing is considered. In sect. 4 we introduce observables which are 
sensitive to sin^ 6'i3 and sign[Am^3]. In sect. 5 we calculate these observables for normal 
and inverted hierarchy and different intervals of sin^ 6'i3, identifying the regions in the space 
of these observables in which the mass hierarchy and sin^ ^13 can be determined. In sect. 
6 we elaborate on the method of the high energy tails of the spectra produced in the 
detectors by neutrinos and antineutrinos. Sec. 7 is devoted to the study of the possibility to 
restrict sin^ 6*13 using the time variations of signals induced by the shock wave propagation. 
Conclusions are given in sect. 8. 



2 Fluxes, density profile, neutrino mass spectrum 
2.1 Properties of supernova neutrino fluxes 

In this section we summarize our present knowledge of the neutrino fluxes from the collapsing 
stars. 

At a given time t from the core collapse the original flux of the neutrinos of a given 
flavor, can be described by a "pinched" Fermi-Dirac (F-D) spectrum, 

L E'^ 

Fl{E, T„ r7„, L„, D) = ^^j^^^'^p^^^^-^ ^e/t^-^^ + 1 ' (1) 

where D is the distance to the supernova (typically D ~ 10 kpc for a galactic supernova), 
E is the energy of the neutrinos. La is the luminosity of the flavor i/^, and Tq represents 
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the effective temperature of the gas inside the neutrinosphere. Supernova simulations 
provide the indicative values of the average energies |27j : 

(E,) = (14 - 22) MeV, {E,)/{E,) = (1.1 - 1.6), {E,)/{E,) = (0.5 - 0.8), (2) 

and the typical value of the (time-integrated) luminosity in each flavor: ~ (1 — 5) ■ 
Iq52 gj^gg_ fjj^g intervals in eq. Q include the results of recent Monte-Carlo simulations 
|27j . in which a difference between (Er^) and (Eg) of about 10% is given as typical, though 
differences as small as few per cent are not excluded. We will briefly comment on this latter 
case in the discussion of our results. The z/^ and u^- {p^ and Vr) spectra are equal with 
good approximation, and therefore the two species can be treated as a single one, {i^x)- 
Conversely, small differences exist between the energy spectra of the non-electron neutrinos 
and antineutrinos, Vx and due to effects of weak magnetism, [2H]- In particular, the 
difference 

5x = Tx-Tx = 0.4 - 0.5 MeV (3) 

is found to be valid with ~ 20% accuracy over a large variety of physical situations |28j . 
The luminosities of all neutrino species are expected to be approximately equal, within a 
factor of two or so [TJ EOj : 

Le/L, = (0.5-2), L,/Lx = (0.5-2), L, ^ Lx. (4) 

The pinching parameter rja takes the values 

r/e~0-3, r/g~0-3 r/^ = r/, ~ - 2. (5) 

Notice that the and z/g spectra may have stronger pinching than the other flavors [211 EO] • 
In eq. (0) the normalization factor -^3(770) equals: 

poo ^3 

^3(^a) = / ^x-r, , • (6) 

Jq + 1 

In the absence of pinching, rja = 0, one gets -^3(0) = 77r^/120 ~ 5.68. The average energy 
{Ea) of the spectrum depends on both Ta and 77^; for ?7q, = we have {Ea) — 3.15 T^. 

The luminosity, temperature and pinching of the neutrino flux vary with the time t; Time 
dependence can be different for different stars j2ZI- If these variations occur over time scales 
larger than the duration of the burst, r ~ 10 s, the integrated z/q flux is can be described 
by the expression ((H) but with smaller rja (integration leads to widening of spectra). 



2.2 Matter profile 

In our study of neutrino conversions inside the star we use the following matter density 
proflle of the star: 

p(r) = 10^3 c ( ' g . cm-3 , (7) 
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with C~l-15|SIll2niE3l2ni- Forp~(l - 10^) g ■ cm-^ expression © provides a good 
approximation to the calculated matter distribution during at least the first few seconds 
of the neutrino burst emission. For p ^ 1 g ■ cm~^ the exact shape of the profile depends 
on the details of evolution of the star, its chemical composition, rotation, etc.. As it was 
recently pointed out after ~ (2 — 10) seconds from the core collapse and bounce, the 
shock-wave propagating inside the star could reach the regions relevant to neutrino resonant 
conversion and modify the observed neutrino signal. 

For three active neutrinos the difference of the v^. and potentials in matter depends on 
the number density of electrons: rig = Yep/niN, where rriN is the nucleon mass and Ye is the 
number of electrons per nucleon. The transitions occur mainly in the isotopically neutral 
region where Yf, = 1/2 with rather good precision. 

2.3 Neutrino mixing and mass spectra 

We study the effects of neutrino flavor conversion in the star and in the Earth in a three 
neutrino framework with mixing and mass splittings allowed by the atmospheric and solar 
neutrino data. 

The neutrino mass eigenstates, Vi {i = 1,2,3), with masses rui are related to the flavor 
eigenstates, (a = e,/i, r), by the mixing matrix U: = ^^Uaii^i- The standard 
parameterization of this matrix is used here (see e.g. [SI])- The atmospheric neutrino data 
determine 

ml-ml = Am^2 = ±(1.5 - 4) ■ 10" tan^ 623 = 0.48 - 2.1 . (8) 

The two possibilities, Am^g ~ ^^31 > and Am^g ~ ^1^31 < 0, are referred to as normal 
and inverted mass hierarchies/ordering respectively and will be denoted as n.h. and i.h. in 
the text. 

The solar neutrino data and the KamLAND results select the large mixing angle (LMA) 
MSW solution with parameters (see e.g. [Hj, [HT)]-p1]): 

ml-ml = Aml-^ = (4 - 30) • 10"^eV^ tan^ 612 = 0.25 - 0.85 . (9) 

The reactor experiments CHOOZ and Palo Verde give the upper bound on 13-mixing |10l 

E]: 

sin^ ^13 < 0.02 . (10) 

3 Mass hierarchy, Ues and conversion effects 
3.1 Permutation factors 

The fluxes Fe and Fg of the electron neutrinos and antineutrinos in the detector can be 
expressed in terms of the permutation parameters (1 —p) and (1 —p) and the original fluxes 
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as follows [S]: 



Fe = vF! + (1 - v)F'^ = F^ + (1 - - F°) 

F, = pF^ + il-p)F',=F^+il-pm-F^) . 



(11) 
(12) 



The factors p and p are the total Ue and z/g survival probabilities which describe the conver- 
sion effects inside the star, in the intergalactic medium, and, if the Earth is crossed, in the 
matter of the Earth. As can be seen in eqs. (fTTl [T^ . the conversion effects are proportional 
to the difference of the original (z/e) and (z/^.) fluxes. 

As they propagate in the star, the neutrinos undergo two MSW resonances (level cross- 
ings) : 

(i) the first resonance (H) occurs at higher density, p ~ 10^ g • cm~^(10MeV/-E'), it is 
governed by the atmospheric mass squared splitting, Am32, and by the ^13 angle. The 
H-resonance is in the neutrino (antineutrino) channel if the mass hierarchy is normal (in- 
verted). The probability of transition between the eigenstates of the Hamiltonian (jumping 
probability) in this resonance, Ph, can be written for the density profile (jZj) as jS]: 



The probability has been obtained on the basis of the Landau-Zener (LZ) formula 
which holds for the linear density distribution. It can be checked that for small mixing (as 
it is the case here) the formula works well for the density profile (jZj) E! ^ IH]. 

(ii) A second level crossing (L) determined by the "solar" parameters Q, happens at lower 
density, p ~ (30 - 140)(10MeV/£') g ■ cm~^. For Am^^ and 6^12 in the LMA region the 
propagation through this resonance is adiabatic for all possible values of C 

Let us first consider the effects of conversion in the star with no Earth crossing. We will 
use here the approximation of "factorized dynamics" in the resonances, which means that 
the dynamics of level crossing in the two resonances are independent and the total survival 
probability is the product of the survival probabilities in each resonance (see |3] for details). 
In this approximation the survival probabilities p and p equal [Hj: 




(13) 



(14) 



P - PH\Ue2\' + {I 
P^ \Uel\^, 



PH)\Ue3 



(15) 
(16) 



for n.h. and 



P^ \Ue2\ 

p - PnlUeil^ + (l 



PH)\UeZ 



(17) 

(18) 



5 



for i.h.. The transition from normal to inverted hierarchy cases corresponds to the inter- 
change of the permutation factors, p ^ p, and of the indices 1 ^ 2 in the mixing parameters. 
Using the standard parameterization of the mixing matrix, according to which 

|t/el|2 = cos' discos' ^13, | f/e2 T = sin' ^12 COs' ^13, | f/e3 1 ' = sin' ^13 , (19) 

we can rewrite the survival probabilities in the following form 

p ~ P^^ sin' + [1 - Ph{1 + sin' ^12)] sin' ^13 , (20) 
p~cos'ei2(l-sin'ei3) , (21) 



for n.h. and 



for i.h.. 



p ^ sin' ^12(1 -sin' ^13) , (22) 
P^Ph cos' + [1 - Pnil + cos' ^12)] sin' ^13 , (23) 



The permutation factors depend on sin 613 via the jumping probability Ph = Pff (sin ^13) 
and explicitly via the mixing parameters. A further dependence of the permutation factors 
on ^13 is given by small terms which are neglected in our approximation. They arise from: 

1) generalizing the LZ expression for Ph'- a more precise double exponential form of 
the jumping probability Ph (see e.g. |15) gives sin' 6*13 corrections to the Landau-Zener 
form. These corrections appear however with very small coefficients and therefore can be 
neglected. 

2) relaxing the approximation of factorized dynamics: this could give linear (in sin6'i3) 
corrections of the form sin 6'i3Am2]^/Am3j^. For the large ^rn^i part of the LMA region we 
have I Am|;^ / Am';^ I ~ 0.1, which is large enough so that sin^i3|Am';^/Am';^| ~ sin' ^13 and 
the linear terms can not be neglected with respect to the sin' ^13 ones. In contrast, taking 
the best fit values of the mass squared splittings one gets |Am2]^/Am|]^| ~ 10"', therefore 
the sin' 9rs terms dominate for sin' ^13^ 10~^. A detailed study of the deviation from the 
factorization will be given elsewhere 



Two limits are important: 

1) . If 613 is very small so that Ph = 1, both mass hierarchies lead to the same result: 

p ~ sin' 6^12, p — cos' 9i2 ■ (24) 

2) . For large enough 6*13 (in practice for sin' 6*13^ 10~^ see the next subsection) when 
Ph = we have 

p ^ sin' 9i3 , p cos' 6'i2(l — sin' 9rs) ^ cos' 612 (n.h.) , (25) 

and 

p ~ sin' 6^12(1 - sin' 6^13) ^ sin' Ou , p^ sin' ^13 (i.h.) . (26) 

Let us underline that in this limit uncertainties related to the density profile or with fac- 
torized dynamics disappear. The conversion is strongly adiabatic and the results ()25| l26|l 
are exact. Notice also that from the point of view of present oscillation results (large or 
maximal 12 and 23 mixings) such a possibility looks rather plausible. 
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3.2 Ph and sin^^is 



The jumping probability Ph is shown in fig. ^ (upper panel, central lines) as a function 
of sin^ for two extreme neutrino energies and the parameters lAmggI = 3 ■ 10~'^eV^ and 
C = 4. As follows from eq. (fTSj) . a constant value of Ph corresponds to sin^ oc E"^^^. 
This determines the shift of the curve Ph in Fig. Q along the axis sin^ ^13 with change of 
energy. 

From the figure it follows that, for all neutrino energies relevant for observations, the 
conversion in the H resonance is adiabatic {Ph = 0) for sin^6'i3^ 4 ■ 10~^. The adiabaticity 
is maximally violated {Ph = 1) for sin^ ^13^ 4- 10"^. In the intermediate region Ph depends 
significantly on sin 6'i3, decreasing from 1 to as sin ^13 increases. 

Let us consider the uncertainties in the relation between Ph and sin^ On- The uncertainty 
in Ph due to a variation of C by factor (1/4 - 4) is shown in fig. ^as a shadowed band. As 
given in eq. (fT^ . a constant value of Ph corresponds to sin^ ^13 oc C~^/^. Thus, for a given 
Ph, the variation of C corresponds to a variation of sin^ 6*13 by factor 2.6. 

Furthermore, the profile ((Zj) is a simplification and deviations from cubic dependence 
can appear if [IH] : 

1. the energy transfer in the star is not purely radiative but also convective. 

2. the luminosity L is not constant with the radius r ^6] . 

3. the opacity k has a spatial dependence: k = k{r). 

In general, one can fit a realistic profile which appears from numerical models of super- 
nova progenitors by p oc r~"'^^^ It can be seen that n = 3 provides the best power law fit of 
the realistic profile, however local deviations can be significant with variations of index in 
the interval n = 1 5 pRTj. 

As discussed in ref. jS], if the position of the resonance, r^es, is kept fixed, the adiabaticity 
parameter 7 scales as n~^: 

|Amf3| sin^2gi3 1 _ 2\Amlr^\rres sin^ ^13 .^^s 

^ ~ 2E cos 2^13 \{dn/dr)/n\res ^ E n ' ^ ' 

So, the uncertainty n = 1 ^ 5 corresponds to a factor ~ 2 — 3 uncertainty in sin^ 61^ with 
respect to the n = 3 case. 

Thus, the absolute limitation in determination of sin^ 61^ from measurements of Ph, 
can be characterized by the factor 1/3 - 3, unless the knowledge of the star profile will be 
substantially improved. In principle, the uncertainty on the density profile can be reduced 
if the progenitor of the star is identified and detailed observations of the supernova optical 
signal (light curve, etc.) is done. 

The uncertainty due to the error in Amg^^ is numerically less important. We take it to 
be of ~ 20%, in agreement with the expected precision of near future measurements |48j . 
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Ph as a function of energy for different values of sin^ ^13 is given in the bottom panel 
of fig. H In the observable part of the spectrum, E = [5 — 70) MeV, a strong dependence 
of Ph on E, and consequently the strongest distortion of the energy spectrum, is expected 
if sin^6'i3 = (1 — 5) ■ 10~^. For sin^6'i3 = 4 ■ 10~^ the probability Ph increases with E 
by factor of 5 in the observable part of spectrum. Effects due to this strong change could 
give the possibility to probe 6*13 is this region. The energy dependence is very weak for 
sin^ ^13 > 4 ■ 10"^ and sin^ ^13 < 4 ■ lO'^. 

3.3 Three regions 

According to eqs. (j2()|) - (pT| for the normal mass hierarchy, the probability in the antineu- 
trino channel, p, does not depend on the neutrino energy and only very weakly depends on 
^13 via the sin^ ^13 term. In contrast, for the neutrino channel the permutation factor, 1 — p, 
depends on the energy via the Ph parameter and on the ^13 angle both explicitly (the sin^ ^13 
term in eq. (HOI)) and implicitly through Ph (eq. (O)- This dependence is illustrated in 
the upper panel of fig. [21 for different values of the neutrino energy, tan^6'i2= 0.38 and the 
other parameters as in fig. ^ According to the expression (pn|l and to the figure we can 
distinguish three regions: 

(i) Adiabaticity breaking region: 



For these values of sin^ ^13 one gets Ph — 1 and the first term in eq. (j2(H) dominates. For 
sin^ 6'i3^ 0: p ~ sin^ 9i2 with very good approximation, independently of the values of 
sin^ ^13 and of the neutrino energy. In this region the conversion in the H resonance has 
little effect and the permutation is due to adiabatic conversion in the L-resonance. 

(ii) Transition region: 



In this region Ph takes intermediate values between 1 and (see fig. Q). Still sin^6'i3 cor- 
rections to the permutation factor are negligible, so that p ~ P//sin^6'i2; the permutation 
(1 — p) increases with sin^ 6*13, following the jumping probability Ph- 

(iii) Adiabatic region: 




(28) 




2/3 



(29) 




(30) 
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Here Ph — and therefore the two terms in eq. ()20|) are comparable near the border of the 
region and with increase of sin^ ^13 the second term dominates. The survival probability 
p has a minimum at sin^ ^i3~ 10~^, corresponding to nearly maximal permutation. For 
sin^6'i3^ 10~^ the term of order sin^ ^13 in the permutation factor starts to dominate and 
p ^sin^ ^13. 

For inverted hierarchy similar results are found with the substitution v ^ v: p has only 
a (small) explicit dependence on sin^6'i3, while p depends on sin^6'i3 both explicitly and 
implicitly (see eq. ((221) )• The same three regions discussed here can be identified according 
to the adiabaticity character of the H resonance. The explicit dependence of p on sin^ ^13 
dominates for sin^ ^13^ 10~^, where p ~sin^ ^13. These features are shown in the lower panel 
of fig. El 

The effects of the explicit dependence of p (or p) on sin^ ^13 in observable signals are 
smaller than ~ 2% and it will be very difficult to study them due to larger experimental and 
theoretical uncertainties. However, their identification may be within the reach of the next 
generation large volume detectors like HyperKamiokande |^ , UNO |^ 1^ and TITAND 

m 

If the accuracy of the experiments is not better than a few % one can neglect the explicit 
dependence of the probabilities on ^13. In this case the permutation factor and therefore 
the observable effects depend on ^13 only via the jumping probability. So one can solve the 
problem in two steps: 

- measure Ph immediately from experiments 

- use the relation between Pjj and ^13 to determine sin^ ^13. 

Summarizing, according to fTTjl , (fT^ and - (j2Sl) the effects of 6*13 consist of 

• change of the degree of permutation; 

• distortion of the energy spectrum due to dependence of the jumping factor on energy. 

For normal hierarchy a change of the neutrino permutation factor, (1 — p), is smaller 
than 30%: from cos^ 6*12 ~ 0.73 at very small 6*13 to ~ 1 at ^13 > 10~^. The change of 
antineutrino factor is negligible. Notice that, for any value of ^13, the conversion inside the 
star leads to the composite spectrum and 6'i3 changes the level of this compositeness. 

Thus, to determine 6'i3 in the case of normal hierarchy one needs to distinguish between 
completely permuted z/g spectrum, with Fg ^ and strongly (3/4 or more) permuted 
spectrum. In the first case one expects a signal which corresponds to the Fermi-Dirac 
spectrum, whereas in the second case, a composite spectrum appears with the dominant 
hard component. 

Clearly, to get an information about ^13 one needs, in general, to know the original 
neutrino fiuxes with better than 30% accuracy (including also substantial uncertainties of 
the detection procedure). 
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For inverted mass hierarchy the change of the antineutrino permutation due to ^13 is 
more significant: (1 — p) increases from sin^ 612 ~ 0.27 for very small ^13 to about 1 for 
^13 > 10"'^. However, in this case the effect of permutation is suppressed due to smaller 
difference of the original Ue and Ux fluxes. The change of neutrino permutation is negligible. 

The permutation of and (or z/g and u^) spectra, and therefore the observed neutrino 
signal, are mostly sensitive to sin^ 6*13 in the transition region; any dependence on sin^ 613, 
outside this interval is negligible. In other words. 



• for sin^0i3> 10~^ the effect of 1-3 mixing is strong in neutrino (antineutrino) chan- 
nel if the mass hierarchy is normal (inverted). This makes it possible to determine 
the hierarchy of the neutrino mass spectrum, while it is difficult to measure sin^^i3. 
Observations of corresponding conversion effects will allow to put a lower bound on 
sin^ 013. 

• in the region sin^ ^i3~ (3 ■ 10"^ — 3 ■ 10^^) measurements of sin^ ^13 are possible, or at 
least both upper and lower bounds on sin^ 6*13 can be obtained; values of sin^ ^13 in this 
region are at least one order of magnitude below the sensitivity of planned terrestrial 
experiments; 

• if sin^ 0i3< 10^^ no effects of sin^ 613 should be seen and the observations are insensitive 
to the mass hierarchy. It follows that in this case the hierarchy can not be probed, 
while an upper bound on sin^ 613 can be obtained. 



4 Energy spectra. Observables 
4.1 Detected signals 

Let us consider the effects of 13-mixing on the energy spectra of the events induced by z/g 
and Ue in the terrestrial detectors. 

The number of charged current (CC) events produced by the z/e-flux with electrons 
having the observed kinetic energy Ee equals 



dE, 



e 



^ ^ J dE'^n{E,, E'JS{E'^) J dEF^{E)^^^^^ , (31) 



where E'^ is the true energy of the electron, A^^- is the number of target particles in the flducial 
volume and S represents the detection efficiency. Here da{E'^, E) / dE'^ is the differential cross 
section of the detection reaction and TZ{Ee, E'^) is the energy resolution function. The Ue 
flux in the detector. Eg, is given in Eqs. ()lipi2|l . 

An expression analogous to (|!TT|l holds for the Ue flux. Eg. 
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The spectrum of the v^. induced events can be measured in water cherenkov detectors, 
hke SuperKamiokande (SK). The relevant CC processes are: 

Vf. +p — ^ 77, + e"*" , (32) 
i/e +0 ^ F + e" , (33) 
i^e +0 -> AT + e+ . (34) 

They are isotropically distributed and essentially indistinguishable from each other. The 
process fl32j) largely dominates the event rate because of the much larger cross section: it 
produces ~ 10^ events for a supernova at D = 10 kpc jHSllHl- Events from z/g scattering on 
electrons can be distinguished, and therefore subtracted, because of their good directional- 
ity; other + e~ scattering processes are neglected due to their small cross section. 

The spectrum of the z/g induced events can be measured in the heavy water cherenkov 
detector SNO. The relevant CC reactions are: 

i/g +d ^ _p + p + e" , (35) 
Ve +(i ^ n + n + e"*" . (36) 

For a typical galactic supernova at D = 10 kpc (and LMA solar neutrino parameters) the 
processes pHjl and give ~ 300 and ~ 150 events respectively |S11IH1- In the near future, 
after the instrumentation of SNO with the ^He neutral current (NC) detectors, the events 
from (jHUj) will be distinguished with ~ 75 — 80% efficiency due to the capture of at least 
one neutron on ^He or on deuterium in coincidence with the detection of the charged lepton 

The light water volume (1.4 kt) of SNO will give signals analogous to those discussed 
for SK (eqs. (jSD-dSl). 

In what follows we consider z/g and z^g- events at SNO and SK, assuming that they can 
be well distinguished and their energy spectra can be separately measured. Our methods of 
analysis can in principle be extended and applied to other types of detectors, e.g. scintillator 
and liquid argon experiments. For a discussion of those in the context of supernova neutrinos, 
we refer to the papers by the LVD jSHl, ICARUS [SIlEHl and LANNDD ^ collaborations. 

4.2 Observables 

Using the differential spectra defined in (jHT|) we can introduce a set of observables which 
are sensitive to the effects induced by the 13 mixing and depend on the mass hierarchy. 

1. The average energy. We define the average energy of events induced by z^e in a 
detector as 
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where 

Ntot= / dE,^ (38) 
Je,^ 'dE, ^ ^ 

is the total number of events above the threshold energy Eth- We take Eth = 5 MeV for 
neutrinos and Eth = 7 MeV for antineutrinos. 

Notice the different energy dependences of the detection cross sections: a oc E^'"^^ for 
z/e in heavy water [^0], while a oc E"^, with negative corrections at high energy, for 1/^ +p 
reaction in water jHll 1^2] (see also the recent calculation and discussion in jHHj)- This 
influences the observables. 

The observed energy spectrum of events includes - in addition to the energy dependence 
of the neutrino flux - the energy dependence of the detection cross section and of the 
efficiency and energy resolution function of the detector. As a first approximation, the 
effect of the efficiency and energy resolution can be neglected and the cross section can be 
taken as (t{E) oc E"^. For simplicity we can also put the threshold of integration to be 
zero. Under these conditions we find the average energy of the observed spectrum of events 
induced by the flux ([T]): 

Li3(-e'?") 

where cr{r]a) is a function of the pinching parameter expressed in terms of the polylogarithmic 
functions: 

OO l. 



k=l 

For unpinched spectrum we find a{0) = 3l7r^/(5670C(5)) ^ 5.07. The quantity a{ria) 
increases with rjc,: e.g. a{2) ^ 5.33 (this can be compared with characteristics of the 
Fermi-Dirac spectrum: a{0) = 3.15 and a(2) ~ 3.6). 

If the neutrino flux in the detector is an admixture of the colder z/g (or u^) and the hotter 

(or Ux) original fluxes, due to conversion effects, its average energy takes an intermediate 
value with respect to the average energies of the two component spectra. It follows that 
conversion effects can be probed by studying average energies. 

The analytical calculation with non-zero threshold gives a more complicated result, which 
is less transparent and is not very different numerically: with Eth = (5 — 7) MeV the 
difference in the average energies is typically 1 - 2 %. 

An important test parameter is the ratio of the average energies of the observed spectra 
of Ue and Ue events: 

- - i • 
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2. The width of the spectrum. The relative width of the spectrum can be characterized 
by the dimensionless parameter 

AE 

r = ^ , (42) 



(E) 



where AE is defined as: 



AE^^{{E-{E)y) = ^{E^)~{Ey . 



(43) 



Here (E) is the average energy and {E"^) is the average of the energy squared. The latter is 
defined by the integral (pTfjl with the substitution Ee E^. 

Measurements of the width will allow to test the compositness of the observed spectrum. 
Performing the integration from Et^ = 0, we find that the relative width of the distribution 
induced in the detector by a Fermi-Dirac neutrino spectrum does not depend on the 
temperature and is determined by the pinching parameter only: 



6Li7(-e''-)Li5(-e''"; 
5 (Li6(-e'''^))2 



.1/2 



- 1 



(44) 



This gives TpulO) ^ 0.44, while r^£)(3) ~ 0.38, corresponding to a ~ 15% narrowing of 
the spectrum with respect to the no-pinching case. To get an idea about effect of the 
compositeness on the width, we give its expression for small energy difference between the 
two original neutrino fluxes: e = {Ex)/{Ee) — 1 <C 1 and for equal pinching, rje = rjx = rj, 
and equal luminosities in the two flavors: 



T{ri) ~ Tpniv) 



(45) 



Moreover, this result is valid if the survival probability, p, is independent of the neutrino 
energy. Eq. (jlKjl shows that r(?7) > rFD{v)y if the permutation is partial (0 < p < 1). In 
presence of strong pinching the width of the composite spectrum could be smaller than the 
width of non-permuted unpinched spectra, i.e. r(?7) ^ Tp£){0). For instance, for 77 = 3, e = 
0.4 and p ^ sin^ 6'i2 ^ 1/4, eq. (j45|l gives r(3) ^ 0.4, to be compared with TpDiO) — 0.44. 
It follows that an observation of F > 0.44 will testify for the composite spectrum and partial 
permutation, while F < 0.44 will testify for pinched spectra without conclusive information 
on the amount of their permutation. 

We introduce also the ratio of the widths of the observed energy spectra of z/g and 
events: 

rr = E . (46) 



3. Numbers of events in the high energy tails. The physics and analysis become simpler 
for energies substantially above the average energy of the spectrum: ii^ > (2 — 3) ■ (E). We 
will refer to these parts of the spectra as the tails. 
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Let us define the number of CC events induced by v^, in tlie detector witli visible energy 
above Ei^: 

r"^ (IN 

N,{E >El)= dE,-^ . (47) 

Similarly, one can define Ng{E > El) as the number of events induced by z/e above the 
energy El- For definiteness in our numerical estimations we will use El = 45 MeV and 
El = 55 MeV. A detailed discussion of the prescription for the choice of these thresholds, 
as well as the dependence of the results on their values, is given in sect. 6.2. 

The following analytical study is useful for the interpretation of results. If El ^ {E), 
the original flux of a given flavor, F^, above the cut is well approximated by the Maxwell - 
Boltzmann distribution: 

^ const ^E^e-^/^" , (48) 

a 

where the normalization constant depends on pinching parameter. 

An estimate of the number of events is given by the convolution of this expression 
with the detection cross-section and the detection efficiency; the energy resolution function 
can be neglected in a first approximation. At high energies the detection efficiency does not 
depend on energy, and therefore factors out of the integration. Taking the cross section as 
a{E) oc E"^, from (jlHI) one gets: 



oo 



^ const LaTa / dxx'^e 

= const LaTaC-^'^^^'-PiEL/Ta) , (49) 

where 

P{x) =x^ + Ax^ + 12x^ + 2Ax + 24 . (50) 
For very high cut, EL/Ta ^ 1, eq. ipUj) has the asymptotic behavior: 

iV° ^ const ^£;le-^^/^" . (51) 

a 

Notice that the result ()49|1 has an indicative characted only, especially in view of the fact 
that the assumption a{E) oc -E^ is a rather crude approximation. In spite of that, however, 
the form (jmj) turns out to be in acceptable agreement with the numerical results, as will be 
discussed later. 

Let us define the ratio of the neutrino and antineutrino events in the tails: 

R,.u{EL,EL)^^^^^. (52) 

taUK L, LJ N_^E>El) ^ ' 

This turns out to be very a powerful test parameter of the conversion induced by the 13- 
mixing. 
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4- Low energy z/g and Vf, events. Similarly we introduce the numbers of events induced 
by Vf. and i>e with visible energies below E'^^ and respectively. In what follows the equal 
values E'l^ = E'^^ = 25 MeV are taken for illustrative purpose, with the low energy thresholds 
Eth = 5 MeV for the neutrino events and Eth = 7 MeV for the antineutrino events. 

We introduce also the ratio of numbers of the low energy neutrino and antineutrino 
events: 

S.^^. (53) 

In what follows we will calculate predictions for these observables depending on the type of 
mass hierarchy and interval of ^13. 



5 Identifying extreme possibilities. Scatter plots 

One can perform the analysis of the supernova data in two steps: 

1) . Resolve the ambiguities related to hierarchy "normal - inverted" and to value of ^13: 
"large - small" . At this point the bounds on ^13 can be obtained only. 

2) . Once the hierarchy is determined and bounds on ^13 are found, one can proceed with a 
detailed analysis of data to measure ^13. 



5.1 Extreme cases 

In this section we consider the first step. We will update the analysis of ref. [S] taking the 
LMA solution with the most plausible values of the oscillation parameters. As in [S], we 
will denote by large ^13 the values of ^13 which lead to Ph ~ 0. The interval is determined 
by eq. (jHUjl . We refer to small 6*13 to indicate values which satisfy eq. (j2HI) so that Ph ~ 1- 
Then there are three extreme possibilities: 



• A. Normal hierarchy - large 6*13: In this case the permutation factors for neutrinos 
and antineutrinos equal: (1 — p) ~ 1 and {1 — p) ~ sin^ 6^12 ~ 1/4 (see eq. (j^ ). One 
should observe unmixed completely permuted (and therefore hard) z/g-spectrum, and 
composite weakly mixed (~ 1/4) i^g-spectrum. 

These features can be quantified in terms of the avarage energies and widths of the z/g 
and z>e fluxes in the detectors: {E{ue)), {E{ue)) and T{u). These quantities are defined 
analogously to the corresponding ones for the observed spectra of events, eqs. (jHTjl 
and ()42|) ^ where AE is determined similarly to the width of the observable spectrum 
P3|l . In general, we get: 

(Eju, )) > {E{u,)), r(z/,) < r(z^e). (54) 

^We mark that, in contrast with sec. 14.21 the present discussion refers to the average energies and 
widths of the spectra of the neutrinos arriving at Earth, and not to the spectra of the events induced by 
the neutrinos in the detectors. 
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The inequality of the widths can be violated in the particular case in which compen- 
sations occur between pinching and permutation effects. Indeed, as discussed in sec. 
14. 2[ equal or slightly larger width could be realized if the original spectrum (and 
therefore the spectrum arriving at Earth) has no pinching, while the original z/g 
spectrum is strongly pinched and the difference between the z/g and average ener- 
gies is small (see eq. ()45|) ). We consider this arrangement rather exotic, since similar 
pinching is expected in the different flavors P?] . 

The Earth matter effect is expected in the antineutrino channel and not in the neutrino 
channel (unless significant difference of the and z/^ fluxes exists jHl] ) . 

If this case is identified we will be able to determine the mass hierarchy and put lower 
bound on sin^ Q\j,. 

• B. Inverted hierarchy - large ^13: The permutation factors equal (1 — j>) ~ cos^^i2 ~ 
3/4 and 1 — p ^ 1. In this case Ue will have unmixed completely permuted F-D spec- 
trum, whereas the z^g-spectrum, should be composite with rather strong permutation. 
Consequently, 

(EiUe)) < {E{Ve)). r(z/e) ^ r(z/e). (55) 

The Earth matter effect is expected in the neutrino channel. 

If this is realized, we will conclude on the mass hierarchy and put a lower bound on 
sin^ Oi^. 

• C. Normal hierarchy - small ^13. Inverted hierarchy - small 9i^. 

These two cases lead to identical consequences: composite z/g-spectrum with cos^ 9i2 ~ 
3/4 mixing (permutation), and composite z/g-spectrum with small sin^ ^12 ~ 1/4 mix- 
ing (permutation). 

Since the permutation is stronger in the neutrino channel one expects: 

{E{v,))>{E{v,)), r(z/e)>r(z7,); (56) 

however these inequalities are not strict since the permutation effects are partially 
compensated by the fact that the original z/g spectrum is softer than the z/g spectrum. 

The Earth effect is expected both in the neutrino and antineutrino channels. 

In these case one can put an upper bound on sin^ ^13 only and the hierarchy will not 
be identified. 

Studies of the z^e-signal only allow, in principle, to disentangle the case A, for which the 
spectrum is hard and of the Fermi-Dirac type, and the cases B, C which lead to the same 
composite spectrum with permutation cos^ 6'i2 ~ 3/4. 

If the i>e-signal is studied only, one can disentangle the case B, characterized by a hard 
Fermi-Dirac spectrum, from the cases A, C which give the same composite spectra with 
mixture sin^ ^12 ~ 1/4. 
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The comparison of the properties of the v^- and z/g-spectra will allow to distinguish three 
possibilities: A, B, and C. 

As we have marked in sect. I4.2| the composite spectra should be wider than the Fermi- 
Dirac spectrum (unless the parameters of neutrino radiation - luminosity, temperature, 
pinching - strongly change with time during the burst). 

To disentangle the possibilities A-C, one can: 

- search for deviations of the spectral shapes from the Fermi-Dirac form, 

- compare the average energies and the widths of the spectra for neutrinos and antineu- 
trinos, 

- search for the Earth matter effects in neutrino and antineutrino channels. 



5.2 Scatter plots 

The criteria described in the previous section for neutrino spectra become less strict when 

(i) the spectra of the observed events (and not of the neutrinos) are considered (here the 
difference of the energy dependences of the neutrino and antineutrino cross-sections plays 
an important role); 

(ii) uncertainties in the original neutrino fluxes are taken into account. 

(iii) uncertanties in the 12-mixing are included. 

Let us recall that in the cases of very small [Pu = 1) or large [Pu = 0) sin^ ^13 the 
uncertainties on the C parameter and on lAmgj^l have no effect on the physics. 

The number of unknown parameters which describe the original spectra is very large. 
For this reason we have constructed scatter plots of the observables using the following 
procedure: 

1) . We define the space of the parameters over which we perform scanning in the 
following way: the average energies, luminosities and pinching parameters of the original 
neutrino fluxes are taken in the intervals 0, 0, Q and (jSj). We assume that tan^ 6*12 will 
be known with ~ 10% accuracy and, as an example, use the interval: 

tan^ ^12= 0.342 - 0.418. (57) 

2) We take a grid of points in this parameter space. Depending on the case under 
investigation, the spacing of this grid is chosen conveniently, with a corresponding number 
of points between ~ 570 and ~ 10^. The number of points used was smaller for the cases A 
and B. Indeed, the scenario A predicts complete permutation of the fluxes in the neutrino 
channel and, as a consequence, the z/g flux at Earth does not depend on the original flux 

(see eq. flllj) ). It follows that a scanning over the parameters of this flux can be avoided, 
resulting in a smaller number of points. An analogous argument applies for the case B and 
the antineutrino channel: no scanning of the parameters of the original flux is needed, 
since this flux cancels from the calculations (eq. (fT^ ). In contrast, for the scenario C, 
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all the original neutrino and antineutrino fluxes are relevant and a full scanning of all the 
parameters (with a larger number of points) is necessary. 

3) For each point of the grid we calculated the observables r^, r-p, Rtaii, S. The cal- 
culations are performed for z/g events at SNO (from the process (jSSl)) and Ue events at SK 
(from the reaction ()32|l ). We chose the energy cuts E'^ = E'^ = 25 MeV for the calculation 
of S and the thresholds E^ = 45 MeV and E^ = 55 MeV for Rtaii- The prescription for the 
choice of these values is given in sec. 16.21 

4) We do not take into account possible correlations of parameters, scanning the points 
within the intervals independently. This leads to the most conservative conclusions. 

The results are shown in figs. El -El 

The fig. El shows the {Rtot — Rtaii) scatter plot. The ratio of events in the tails, Rtaii, 
is defined in ()52|) and Rtot is the ratio of the total numbers of neutrino and antineutrino 
events: 

Rtot ^ ^ , (58) 

where Ntot for neutrinos is defined in ()38j) . From the figure it follows that in the case B 
the ratio Rtot can not exceed ~ 0.032, while for C Rtot can be nearly twice as large. The 
difference is explained by the different degree of permutation in the various cases, according 
to what discussed in sec. 15.11 In case B the observed antineutrino spectrum is completely 
permuted and therefore maximally hard. This in turn implies the largest rate of z/g events 
due to the larger detection cross section at high energies. In the case C the permutation 
of Ue is partial, giving softer observed Ue spectrum and therefore smaller event rate and 
larger value of Rtot- The i>e spectrum is partially permuted in B and C, while in case A 
it is totally permuted and therefore harder. This implies a larger Rtot with respect to C. 
Numerically the difference between the two cases is small (see figure), because in C the 
amount of permutation for z/g is rather large: 1 — p ~ cos^ 612 — 0.75, making this case close 
to complete permutation. 

The parameter Rtaii has much higher discriminative power than Rtot- Rtaii < 0.06 for 
the case B and practically there is no overlap of the regions A and B. The region C overlaps 
with both A and B, although for Rtaii > 0.22 only A is allowed. 

As follows from the figure, there are certain regions in the Rtaii — Rtot plane where only 
one possibility is realized. For the case of inverted mass hierarchy and large 6*13 (B), the 
region is defined as Rtaii = (0.040 — 0.055), Rtot = (0.02 — 0.03). For the case of very small 
13-mixing (C) there is a band around Rtot — (0.3 ■ Rtaii) with Rtaii < (0.05 — 0.14). The 
normal mass hierarchy case A is the unique possibility for Rtaii > (0.14 — 0.22) where the 
border depends on value of Rtot- 

Clearly it will be easy to identify or discriminate the case B. It might be more difficult 
to disentangle A and C since a significant overlap exists. The overlapping areas correspond 
to similar original fluxes in the different neutrino flavors: conversion effects are smaller for 
smaller difference of the original fluxes, making it difficult to disentangle different scenarios 
of mass hierarchy and 13-mixing. 
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parameter \ scenario 


A 


B 


C 


r^; (events) 


in 17 
i.U - i. ( 


n S7 1 1 
U.o ( - i. i 


no 1 
u.y - i.D 


/ £;l^ilt^U.Lllll(JbJ 


\ 1 


^ 1 


> 1 


' r I " VcliLo I 


U.OU ±.1 


1 Dfi 19^ 


n Q 19^ 


rr (neutrinos) 


< 1 


> 1 


> 1 


Rtail 


0.04 - 0.26 


0.04 - 0.056 


0.035 - 0.21 



Table 1: The ranges of parameters of the spectra of observable events ("events") for three 
possible cases of mass hierarchy and 1-3 mixing. For comparison we present also the cor- 
responding characteristics of the neutrino spectra ( "neutrinos" ) , which are not immediate 
observables. 



The fig. m shows the S — RtaU scatter plot, where the ratio of the numbers of low energy 
events, is defined in fl53|) . This figure is rather similar to fig. El although the 5* parameter 
is more complementary to Rtaii than Rtoti and experimentally 5* and RtaU are independent. 
The case B has higher S and smaller Rtaiu inversely, the scenario A predicts higher Rtaii 
and smaller S. The case C is intermediate between the two. As a consequence, the overlap 
of regions is smaller than in fig. El 

We find the following regions where only one possibility is realized: The case A is unique 
for Rtail > 0.165, the case B is unique for S > 0.017 and Rtaii = (0.040 — 0.055), and C is 
unique in the band S ~ (0.2 — 0.3) ■ Rtaii- 

These features are explained in terms of total or partial permutation, similarly to what 
we have dicussed for fig. El 

The figure El shows scatter plots in the space of the variables Rtaii, fHTjl . and rr 
fl42|) . The corresponding intervals of values of these observables for each of the scenarios 
A, B, C are summarized in the Table H For comparison, the Table displays also the 
expected intervals of the same variables for the neutrino spectra discussed in sec. 15.11 
These parameters are not observable, however they allow to understand the features of the 
observable spectra. 

From the figure and the Table it follows that for the case B the allowed region of pa- 
rameters is relatively small, while for A and C the allowed regions are larger and expand 
over te — 1 — 1.5. Furthermore, rr is mostly larger than 1 in C and mostly smaller than 
1 in A. These results largely follow the expectations for the neutrino spectra, and can be 
interpreted in terms of the amount of permutation according to the discussion in sec. 15.11 

One can see two slight deviations with respect to the predictions for neutrino spectra: 
(1) a significant number of points with r^; > 1 in the case B and (2) an appreciable amount 
of points with rp > 1 for the case A. The latter is in agreement with the possibility that 
the broadening of the spectrum due to compositeness is overcompensated by the effect of 
pinching, as commented in sec. 15.11 Moreover, the following reasons contribute to explain 
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the deviations: 

(i) the small contribution of the original v^. flux in the detected z/g signal. This contri- 
bution is - depending on the value of ^13 - not larger than sin^ 6*12 — 0.25, therefore the 
average energy of the observed v^. spectrum is dominated by the harder component due to 
the original flux. 

(ii) the different energy dependence of the cross sections of the z/g +p and v^. -^d detection 
reactions: the flrst grows like E'^ with suppressing corrections at high energies [^21121; while 
the second has a stronger growth with energy, being proportional to _e2-25 follows 
that in presence of equal and z/g energy spectra, this difference leads to higher average 
energy and larger width of the observed z/g spectrum. 

Large regions of the parameter space exist where only one among the scenarios A, B or 
C is possible. Also regions appear when two of these scenarios are realized. If these regions 
are selected by the experiments, the third possibility will be excluded. From the flgs. EHSl 
we conclude that: 

• The case A is excluded if observations give Rtaii ^ 0.06 and/or rr ^ 1.1. If, in 
contrast, the experiments give Rtaii > 0.22 the normal hierarchy would be identifled. 
This result would be further supported if rp ^ 0.95 is also found. 

The possibility of getting information on Ph and consequently on 6*13 depends on the 
speciflc value of Rtaii (as illustrated in sec. 

• The case B is excluded by large values of Rtaii and large values of the ratio of average 
energies: Rtaii ~ 0.06 and ^ 1.2. The identiflcation of this scenario (and therefore 
of the inverted mass hierarchy) appears difficult due to the almost complete overlap 
with the regions of the case C. Another indication of this scenario would be the result 
ve < 0.95. Again, conclusions on 6^13 depend on the speciflc value of Rtaii- 

• The case C can not be easily identifled. An indication of this possibility would be 
the result te ^ 1.2 and rp ^ 1.05. This combination would exclude B and A, and 
therefore indicate Ph > 0, corresponding to small values of ^13: sin^^i3^ few ■ 10"^ 
(see flg. H}. 

The scenarios in which < Ph < 1 are not shown in the flgures. For n.h. and < Ph < 1 
we expect the allowed region to be intermediate between the regions found for A and C. 
Similarly, for i.h. and < Ph < 1 the region of possible values of parameters is intermediate 
between the regions of cases B and C. For this reason, the conclusions we derived from the 
flgures have essentially an exclusion character and not the character of establishing one 
of the scenarios A, B, C. 

It is clear that the potential of the method we have discussed depends on the statistics 
and therefore on the distance from the supernova. Some estimations are presented in sec. 
16.61 and flg. [HI For a relatively close star {D ^ 4 kpc) the error bars are substantially 
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smaller than the field of points so that the discrimination of different possibilities becomes 
possible. 



It is easy to understand the effect on the scatter plots of choosing more conservative 
intervals of the parameters of the original neutrino fluxes. In particular, for smaller differ- 
ences of the average energies in the different flavors the spectral distortions due to conversion 
become smaller (see e.g. The results approach those expected in absence of conversion 
and therefore are the same in the three secenarios, A, B, C. The corresponding points in the 
scatter plots would be located where the regions for the three cases are closer or overlap. 
Clearly, any sensitivity to the neutrino mixings and mass hierarchy is lost in this situation. 

6 The method of the high energy tails 
6.1 Rtaii and r. 

The uncertainties related to the original neutrino fluxes can be substantially reduced if 

• ratios of the electron neutrino and antineutrino signals are considered; 

• the high energy tails of spectra are used. 

The key point is that in the high energy tails the fluxes of non-electron neutrinos domi- 
nate, and moreover, these fluxes are nearly equal. 

In what follows we study the possibility to use the ratio Rtaii{EL, El) introduced in (j52|l 
to establish the mass hierarchy and to probe 613. As we have seen in the previous sections 
the dependences of the neutrino and antineutrino signals on 6*13 and on the sign of Am^g 
are different. 

Let us consider the predictions for Rtaii in details. This ratio can be written in terms of 
the original neutrino fluxes and the survival probabilities p and p as: 

E,) . ;-<^i+°if^i|P' Q(E,, E,) . (59) 

1 - {p} + a{EL){p} 

Here the brackets () denote the averaging over the corresponding energy intervals (we have 
taken into account the weak energy dependence of p and p). 
The quantity Q is defined as: 

- , _ iV°(E > El) 
QiEL,EL) = , (60) 

where A'^, A^^' are the numbers of events calculated according to eqs. ()47|) and (jHT|) with 
the fiuxes F° and F° respectively. Due to the near equality of the fiuxes F° and the 
astrophysical uncertainties in Q are almost cancelled for an optimized choice of the cuts El 
and El (see sect. Ifj.2j) . 
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In fj59p a and a are the parameters which describe the relative contributions of the 
original v^. and v^. fluxes to the numbers of events above the energy cuts: 

"<^'> " Nm • " WET) ■ 



Let us introduce the ratio: 

which should not depend substantially on the features of the detectors, and is known once 
Rtaii is measured. According to (fH^jl we have: 

1- {p){l-a{EL)) 

The ratio Rtaii is a measurable quantity, and, as we will argue in the next section, Q{El, El) 
can be reliably predicted. Therefore, the study of Rtaii can be reduced to that of the quantity 
r(^L,^L). 



6.2 The factor Q 

The factor Q{El,El) depends on (i) the energy cuts El and El, (ii) the fluxes F° and 
and (iii) on the features of the detection method: cross sections, efficiencies, energy 
resolutions, etc. (for the latter quantities we follow ref. [HS])- Since F° and F° have almost 
equal luminosities, temperatures and pinching parameters, the astrophysical uncertainties 
affect Q only weakly. Moreover, the effect of uncertainties can be further suppressed by 
choosing different thresholds for neutrino and antineutrino events, such to "compensate" 
the difference in the temperatures. 

This can be seen from the approximate form of the ratio Q (eq. ()60|)): 

Q ^ const Ii P(^_l/T'-K (-e,/tmt.) (64) 
T,P{El/T,) ^ ' 

^ const (^^y(^|^ye(-^-/^^+^-/^-) , (65) 

which can be derived from the analytical formulas IKT|l . In eq. ()65|1 we took = 
and ?7a, = rj^. For simplicity we considered the same energy dependence of the and Ue 
cross-sections, (t(z/) oc (t{u) oc E'^, and neglected detector parameters like energy resolution, 
efficiency, etc.. From eq. ()65|1 we see that if = the ratio Q reduces to a constant 
provided that El = El is taken. In the presence of a difference 6x between the and Ux 



22 



temperatures, eq. Q, the dependence of Q on does not cancel exactly. However, one 
can find values of the energy cuts, ^ E^, for which this dependence becomes very weak 
in the relevant range of T^: 

^ ^ for = (5 8) MeV . (66) 

In the approximation dx/T^ -C 1 and El/El — 1 <^ 1, Eqs. and ()66p lead to 

I^-1 + I^f2-3|^V (67) 



El \ El^ 

Taking El = 45 MeV, = 7 MeV and 6^ = 0.5 MeV, from ^ we get El ^ 50 MeV. 
A more detailed numerical study, which takes into account also the different energy depen- 
dences of the cross-sections, give El — 55 MeV. 

The dependence of the factor Q on the temperature is shown in fig. IHl (region between 
the lines), for El = 45 MeV, El = 55 MeV and 4 = 0.35 0.5 MeV. As it appears from 
the figure, over the interval = 5-7-9 MeV, Q varies by about ~ 10% of its value, therefore 
it can be taken as a constant with ~ 10% associated error: 

g = 0.0635(1 ±0.13) . (68) 

As it can be understood from eqs. (j65|l and (|67jl . the error on Q is dominated by the 
uncertainty in the value of 6x'- fixing 6x = 0.5 MeV we get Q = 0.0566(1 ± 0.05) in the 
relevant interval of values of T^.. 



6.3 a and a. Threshold Energies 

The ratios a and a, eq. (jUT|) . depend (i) on the features of the original electron and non 
electron neutrino fluxes, (ii) on the energy cuts El, El and (iii) on the characteristics of the 
detectors consider. Using the Maxwell-Boltzmann approximation, eq. ()48|) . and the results 
()49p50j) . one gets the following approximate expression for a: 

where the dependence cr(z/) oc E"^ has been considered and other detection parameters 
(efficiency, energy resolution, etc.) neglected. The polynome P{x) is given in eq. (jSDI)- 
The dependence of ratio P{El/T(.) / P{El/Tx) on El is weak and cancels in the asymptotics 
El ^Tf,,Tx. In this limit we have: 

From eqs. (j^HI I7n|l it follows that, for a given Tj, the ratio a decreases with the decrease 
of Tg and with the increase of the cut El- We have a — > in the limit El/T^ ^ oo. In 



23 



particular, this implies that the contribution of a to the ratio r (eq. (jHSl)) can be reduced 
to a small correction provided that a high enough cut is chosen. 

Results analogous to eqs. (j69t I70|l and similar considerations hold for a. As a con- 
sequence of inequalities of the average energies, eq. ((21), and of the nearly equality of 
luminosities we have 

a{EL) > a{EL) , (71) 
provided that El and -El do not differ strongly. 

Though approximate, the expression ()69|) is in acceptable agreement with more accurate 
numerical calculations. Scanning the intervals of parameters of the original neutrino spectra 
(I21-I1I) we find the following ranges of a and a: 

a(45MeV) ^ - 0.42 , a(55MeV) ~ - 0.95 . (72) 

As an example, taking the "traditional" scenario with equal luminosities in the different 
flavors and hierarchical (unpinched) spectra - with = 7 MeV, Tg = 5 MeV and Tg = 3.5 
MeV - we get 

a(45MeV) - 0.05, a(55MeV) ~ 0.16. (73) 

In fig. [7|we show the dependence of a and a on the cut energies, E^ and E^, for = 7 
MeV, Le = Lg = Lx, rje = rje = rjx = and different values of the and z/g temperature. It 
can be seen that the decrease of a and a with the cuts is indeed exponential, in agreement 
with eq. However, if the hierarchy of the spectra is not strong (solid lines in the 

figure), for our choice of the cuts a, and especially d, are not negligible with respect to the 
survival probabilities: a ^ a ^ p,p ^ 1. In this case, the effects of the a and a terms 
in r (and therefore in Rtaii) can be reduced by adopting higher energy cuts. According to 
fig. la requiring a ^ a ^ 0.1 implies El ^ 70 MeV and El ^ 100 - 120 MeV (we have 
considered the possibility, not shown in the figure, to have Le/L^ — 2 and/or Lg/L^ — 2, 
corresponding to twice as large values of a and a with respect to what shown in the figure). 
It is clear, however, that the numbers of events in the tails above these energies is strongly 
suppressed, so that no precise measurements of r are possible. 



6.4 Large threshold energy Umit 

For El ^ (Ee), (Eg) parameters a, a are negligible and expression ()63|) becomes 

^ 1 - (p) ' 



(74) 



i.e. the astrophysical and oscillation parts of Rtaii are factorized. 

Inserting the explicit expressions for the permutation factors, eqs. ()20l - and ne- 
glecting ~sin^ 6*13 terms, we find 

r = -4^ - {Ph ) (n.h.), (75) 
sm 0i2 
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where {Ph ) is the value of the jumping probabihty averaged over the energy interval under 
consideration. The ratio r is shown in fig. |H1 as a function of {Ph ) for different values of 
tan2 6'i2. In the limit of {Ph) 1, i-e. very small 6*13, we get 

tan"^ 9i2 

for both hierarchies. In general one finds 

'>l/tan2^i2 (n.h.) 
<l/tan2^i2 (i.h.) ' 

and in the limit {Ph) (large 6*13) 



(78) 




For both the hierarchies r decreases with the increase of the mixing 612- 

Once (i) Rtaii is measured in experiments, (ii) Q is calculated and (iii) tan^ 612 is known 
from the solar neutrino and the KamLAND experiments, the inequalities ()78p can be used 
to establish the mass hierarchy. Assuming that tan^ 6*12 will be measured with 20% accuracy, 
from the figure |H1 (see dash-dotted lines; the central value tan^ 6*12= 0.38 has been taken as 
an example) we find that: 

1) . If r < 2 the inverted hierarchy will be selected. In particular, if r ~ 0.7 even a poor 
accuracy in measurements of r, say 30 — 50%, will be enough. Moreover, the (energy- 
averaged) jump probability will be restricted: {Ph) < 0.3, leading to a lower bound on 
sin^ 6^13. 

2) . If r = (2.0 — 3.5) both types of hierarchy are possible. For inverted hierarchy one can 
put a rather strong bound on the jumping probability: {Ph) > 0.8, which would imply an 
upper bound on tan^ ^^3. In the case of normal hierarchy no bound on Ph appears. 

3) . For r > 3.5 the normal hierarchy with {Ph) < 1 is selected. 



According to eqs. (f7H|) - (ffH|) we have: 



1 



iPH) = -^-r (n.h.), (80) 



sin' 9 



12 



cos^ Ui2 r 

The result for {Ph) can be then transferred into a result for tan^ ^13 using fig. ^ and the 
expressions ((131 EH)- As follows from the fig. [H even for very precise measurements of {Ph) 
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the uncertainty on the density profile will lead to a factor of 3 uncertainty on sin^6'i3. To 
have a sensitivity to sin^ ^13 in the adiabatic region one needs to measure the permutation 
effects at the level 2% or smaller. This looks practically impossible already in view of un- 
certainties on the determination of tan^ 612. 



6.5 The general case 

To have significant statistics in the tails the energy cuts should not be too large. In this 
case a, a can not be neglected and for the ratio r we should use the complete expression 
fj63|) . It can be rewritten in terms of the average jumping probability (Ph) and 6*12 as 

' = tan^^,, + a ^'^^ 

' + (i.h.). (83) 



l + tan2^i2- (PH)(l-a) 
For very small ^13 {Ph = 1) both hierarchies give the same result: 



1 + tan^ 6*120 

'^i = T7\ • (84) 

tan^ 612 + a ^ ' 



Comparing ri with the general expressions ()821 l8o|) we get the inequalities 

r > ri (n.h.), r < ri (i.h.) , (85) 
which provide a test of the type of mass hierarchy. 

In the adiabatic case {Ph ~ 0): 



1 + tan^ 012 



tan 9i2 + a 



(n.h.), (86) 



(j) 1 + a tan^ 6*12 
1 + tan^ 6*12 

These quantities turn out to be the upper (for n.h.) and the lower (for i.h.) bounds on r 



Ti < r < rl^^ (n.h.). 



rj^<r<ri (i.h.). (89) 
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The averaged probability, {Ph)^ can be expressed in terms of measurable quantities as: 

{Ph) = ^ J l + cot^ - r(l + a cot^ ^12)] (n.h.) (90) 
(1 - a) 

{Ph) = _^ [l+tan'gi2- ^(l + atan'gi2)] (i.h.). (91) 



We can define the energy Eav in such a way that 

{PH)=PH{Ea.). (92) 

Using the Landau-Zener formula for Ph we find expression for sin^ 613 in terms of observ- 
ables: 

sin^^,3--2.05.10-^^^j (ikij) ^" = "'^)' 

where (-Pff) are given in Eqs. (fHUI PT| . 

6.6 Results 

Let us first estimate the influence of a and a on the ratio r = r{a, a). We denote here and 
later by and x~ the upper and lower edges of the uncertainty interval of the variable x 
{x = a, a, ...). Using the general expressions for r, eq. (jH2I), and taking into account the 
restriction (f7T|) we find the upper (r^aa;) and the lower (r^m) bounds on r for given values 
of 6*12 and Ph- In the case of the n.h. we get 

X 1 + tan^ 612 - {Ph) tan^ 6^ 

rminiPn) ~ r(0, a^) = — ^ — — , (94) 

tan"^ 9i2 + a+ 

TmaAPH) ~ r(0, 0) = - (Ph) (95) 

sm Ui2 

which coincides with expression in the high energy limit. Here the minimum values a~ and 
a~ have been set to zero for simplicity. 

For inverted hierarchy we find the lower limit: 

1 

l + tan^ eu~{PH){l-a+] 
while, due to the restriction ()7ip . the upper limit has a more complicated dependence 



TrmniPH) = r(0, «+) = ^ _ , , ^^^^7^ ^TT , (96) 



fmax{PH) 



12 



'r(a+,a+) = [1 + a+tan2^i2] [1 + tan2^i2 - (Ph)(1 -a+)] ^ Ph <ta.n^ B 
r(0, 0) = [1 + tan2 6^2 - {Ph)T^ Ph >tan2 6^2 

(97) 



27 



These limits are shown in fig. |H1 for tan^^i2 = 0.38 and = = 1, together with 
predictions for a = 0.17 and a = 0.34 as an example. As follows from the figure, with the 
present knowledge of the original spectra the mass hierarchy can be identified from the tail 
method only if 

r > 2.6 (98) 

will be found. In this case also the upper bound on (Ph), and consequently a lower bound 
on sin^ 6'i3 can be obtained: 

(Ph) < - (99) 

sm 6'i2 

For r < 2.6 neither the hierarchy nor sin^ 6'i3 can be found. For the most plausible 
scenario (normal mass hierarchy, large sin^ 6*13 {P^ = 0) we get r ~ 2, that is, below the 
identification interval. 

More can be said if the hierarchy is identified by some other method {e.g., from the 
Earth matter effect or the shock wave effects, see sees. 16.71 and ITj). Thus in the case of 
normal hierarchy, in addition to the upper bound (jHH), one can put a lower bound on (Ph) 
if r = 0.7-1: 

(Ph) > -r(l + -4^1 ' (100) 

sm^ 6^12 V tan^ 612 J 

and, correspondingly, one gets an upper bound on sin^6'i3. 

If the inverted hierarchy is identified and r = 1 — 2.6 is found, one can put a lower bound 
on Ph (upper bound on sin^ 613) : 

(Pi/) > 1 + tan^ ^12 - 1/r. (101) 

The discrimination power increases with the increase of the low energy cut (i.e., decreases 
with and a~^), however the statistics decreases fast correspondingly, so that such a 
possibility can be realized only in the case of supernova at small distances. 

The features of r are reproduced also in the observable Rtaii, eq. (j^ . We find the 
minimal and maximal values of Rtau = RtaiiiEav, {\Aml2\),C,a,a,9i2,9i3,Q) for a given 
6'i3 taking into account also uncertainties on Am23, (^12, C, Eav For the normal mass 
hierarchy, eqs. (|94p95j) and (jl3pi4j) give: 

Rma.{6l^) = R{E-^, C+, 0, 0, ^13, Q+) 

RrmniOis) = RiE^,„{\Aml2\)-,C-,0,a'-,et2,e,3,Q-) . (102) 
Similarly, from eqs. ()96p97|l we get for inverted mass hierarchy: 

'R{E+,{\Aml2\)-,C-,a+,a+,9^2,0rs,Q^) Ph <ta.n^ 9u 
R{E+,{\Aml,\y,C-,0 ,0 ,er2,^i3,Q+) Ph >tan^ 0^ 

Rrmn{9n) = R{E;^, {\Aml2\)+ , C+, 0, a+, 9+2, ^13, Q"). (103) 



R„ 



'13 j 



R max and Rmin &s fuuctious of siu ^13 for n.h. and i.h. are shown in fig. 1101 In our 
calculations we have used the following uncertainties: 
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• The effective energy Eav practically coincides with the value energy cut, for n.h. 
or El for i.h.. A 10% uncertainty on its value was adopted to be conservative. 

• A 20% error was taken on |Am|2|- 

• The parameter C was taken to be in the interval C = 1 — 15, so that C~^^^ = 0.4 — 1. 

• tan^ 6i2 was assumed to be known with 10% accuracy. 

• the parameter Q was taken according to eq. (jUHj) . 

• a and a were taken in the intervals (ff^ . 

We mark that these assumptions on the uncertainties are consistent with the intervals 
()2l4p ■ which have been taken to generate the scatter plots EE 

Notice that for i.h. Rmin is independent of 6*13. This is due to the fact that Rmm is 
realized in our analysis at a = a"'" ^ 1 (equal permuted fluxes), for which no conversion 
effect appears and the dependence of r on Ph cancels (see eq. (jHBI))- 

To estimate the sensitivity of the method we have simulated three possible experimental 
results for two different distances to the supernova, D = 4 and D = 8.5 kpc and the 
integrated luminosity Lp = 5 ■ 10^^ ergs in each of the six neutrino (and antineutrino) 
species. The error bars correspond to 99% C.L.. 

In absence of Earth crossing, The following relevant intervals for Rtaii are found (see fig. 
Cni): if 

Rtaii ^ Q^/ tan^ 9^^ - Q+/ sin^ (104) 

the n.h. is established and a lower bound on sin^ ^13 is put. If this condition is not fulfilled 
it is not possible to establish the mass hierarchy and conditional bounds on sin^ 6'i3 can be 
obtained. In particular, for 

• R^aii — / tan^ '■ an upper bound on sin^ 6*13 is found in the hypothesis of 
inverted mass hierarchy. 

• Rtaii — Q^'- no information is obtained on sin^6'i3. 

• Rtaii — Q~ cos^ 9^2 -^ Q~- an upper bound on sin^ 6*13 is established in the hypothesis 
of normal mass hierarchy. 

No measurement of 6*13 (that is, no lower and upper bound) is possible. The sensitivity 
to the n.h. (i.h.) increases (decreases) with the increase of 9i2- 

Clearly, conclusions depends on the statistical error on Rtaii as given by the experiments. 
This in turn is determined by the distance to the supernova, by the neutrino luminosities, 
the volumes of the detectors, etc. (see fig. ITn|l . 
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6.7 Including the Earth matter effects 



If the neutrino burst crosses the Earth before detection the regeneration effects in the matter 
of the Earth should be taken into account. The effects increase with neutrino energy and can 
reach 30 - 50% at E > 50 MeV EEl ESI EZI • With these effects the survival probabilities 
p and p become: 

P^PHP2e, P^Ple, (n.h.) (105) 



p~P2e, p^PnPie. (i.h.), (106) 

where P2e and Pie are the probabilities of 1^2 -^Ve and vi conversion in the Earth 

respectively. They have an oscillatory dependence on the neutrino energy and can be written 
in terms of the regeneration factors jS] : 

P2e = Sin^ 6^12 + /r-eg, Pie = COS^ 6'i2 + freg ■ (107) 

According to ()105[ I106[ I107jl the generalization of results to the case of the Earth matter 
effect is straightforward: In the formulas of sect. 16.51 one should substitute sin^6'i2 —>■ 
sin^ 9i2 + freg, cos^ Ou cos^ 6^12 + freg and use the averaged regeneration factor ^: 

_ 1^;;° J dEin{E,,Ei)S{E',) J dEfreg{E)F^{E){da{E',,E)/dE',) 

t = Tnz^ 5 (lOo) 

/;;° / dE'MEe, E'J£{E'J J dEiE)FOiE)ida{E',, E)/dE'J 

for neutrino channel and / (with an analogous definition) for the antineutrino channel. 
From eqs. ()105|) - ()108|1 one finds the following generalization of eqs. ()82p83|) : 



l-(Pg)(sin^gi2 + /)(l-a) 
l-(cos2^i2 + /)(l-a) 

l-(sin2^i2 + /)(l-a) 



(n.h.) (109) 



(i.h.) (110) 



l-(P^^)(cos2ei2 + /)(l-a) 

The quantities / and / depend on Amgi, 9i2, on the direction to the supernova, on the 
Earth density profile, on and and on the energy cuts E^ and El. We calculated the 
values of the averaged regeneration factors using a realistic density profile of the Earth [BH] 
and taking the nadir angles 9^ = 84.3° at SNO and 6'„ = 24.5° at SK [12], with El = 45 
MeV, El = 55 MeV, Am^^ = 5 ■ 10"^ eV^ and tan^ 6112= 0.38. To estimate the uncertainties 
on / and / a 10% error was considered on Arri^i and tan^ 6*12, while (T^) was allowed to 
vary in the - rather large - interval 4.5 -r- 8 MeV ( 5 -^ 8.5 MeV). We get 

/ = 0.21 ^0.24, 7=0^0.06. (Ill) 



^Rigorously, / depends also on a, and not only on the i^x original flux, as given eq. (|108|l . However, it 
can be checked that the corrections due to a are always negligible for the energy cuts in consideration. An 
analogous conclusion holds for the a corrections to /. 
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For different nadir angles, comparable uncertainties on / and / are found. 

Similarly to what discussed in sees. 16.51 and 16.61 using ()110|) we calculate the interval 
of possible values of Rtaii, Rmin Rmaxi taking into account all the uncertainties, including 
those on / and /. The results are shown in fig. ^2 Notice that with respect to the 
no-crossing case, a stronger dependence of Rtaii on 6*13 is seen for n.h., implying a larger 
sensitivity of the data to the normal mass hierarchy. This is related to the fact that the 
Earth regeneration effect is stronger in the v^. than in the z/g channel. 

Let us mark that the observation of oscillatory distortions due to Earth matter effects in 
the energy spectrum of the Ve (t'e) signal can establish the normal (inverted) mass hierarchy 
[3 EH] Once the information on the hierarchy is known, the method discussed here will 
provide a measurement or bound on ^^13. 

6.8 Refining the method 

As we saw in the previous sections appearance of the a-terms substantially reduces the 
identification power of the method (cf. fig 1^ and fTUI) . Let us consider the possibilities to 
refine it. 

1. Further studies of the physics of supernovae can lead to more precise predictions 
for the original neutrino spectra. Furthermore, one can use results of studies of the whole 
neutrino and antineutrino spectra to reduce the uncertainty intervals of the average energies 
and fluxes, which then can be used in the tail method. 

2. One can take larger cut energies and El. With the increase of Ei and E^ the 
a-terms decrease and the picture shown in fig. El will approach that in fig. |S1 For higher E^ 
and El the hierarchy can be identified in a larger interval of r (and Rtau) and the 13-mixing 
can be measured (especially in the case of inverted hierarchy). This will be possible in the 
case the distance to the supernova is relatively small, so that the statistics is high enough. 

3. One can determine a, a or other relevant parameters of the original neutrino spectra 
immediately from the experimental data. This can be done by studying the dependences of 
the numbers of events Ne^El), Ns^El), as well as their ratio Rtaii, on the threshold energies 
El, El- Indeed, the number of events N^^El) depends on El via the quantities and a: 

N,{El) = (1 - p)iV°(i?L)(l -p + pa{EL)) . (112) 

These two contributions in principle can be disentangled, since they are different functions 
of El- As follows from eq. ()7()j) . a has an exponential dependence on El of the form: 
a{EL) = Aexp{—EL/Tex) (with l/T^x = l/Tg — l/T^). A fit of data could provide the 
unknown parameters A and T^^, allowing to fully reconstruct the function a(i?/,). 
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7 Time dependence of the signal: shock-wave effects 



7.1 Shock wave and H- resonance 



As it was pointed out in ref. the shock- wave propagating inside the star may reach 

the region of densities relevant for neutrino conversion during the neutrino emission period 
(10 — 20 seconds). It modifies the density profile of the star, thus affecting the pattern of 
neutrino conversion [221 EHl • 

An indicative description of the time dependence of the matter density profile is given 
in fig. nil (adapted from (HB])- As shown in the figure, at a given time t post-bounce, the 
density distribution presents an overdense region corresponding to the shock- front. The 
outer border of this region (front) consists in a sharp "step" in which the density increases 
from Pi to pf. According to the relative increase of the density in the front, ^, can be 
of about an order of magnitude: 

e = p/M~io. (113) 

Above the shock front the profile coincides with that of the progenitor star. Below the front 
a rarefaction region, a "hot bubble" , is produced. 

According to [SB], the velocity of the shock front, Vg, increases with the distance from 
the center of the star. For distances d = (1 — 5) ■ 10'^ km it changes in the interval 

Vs ~ (0.8 - 1.2) ■ 10^ Km/s, (114) 

and for d > 5 ■ 10^ km, Vg is nearly constant: ^ 1.2 ■ 10^ Km/s. 

The shock wave influences the neutrino conversion when it reaches the resonance layer. 
Given the mass squared splitting Am^, the neutrino level crossing (resonance) is realized at 
the density 

mArAm^ Am^ \ /'lO MeV\ 

Prei 



V2GfE 



cos2»,3^1.4.103g.cm-(^)(H|!X). (US) 



where rriN is the nucleon mass and Gp the Fermi constant and we have taken the electron 
fraction Yg = 1/2. From the fig. p2j) and eq. ((7j) we find the resonance radius ^: 

,„..5.10V.(iy^) (^) , (116) 

and - taking a constant speed Vg with values ()114l) - one gets the time tg ~ dres/vg after 
which the shockwave reaches the resonance: 



■^In ref. the density profile p{d) cx d " is used, with n = 2.4. Here we use n = 3 for simphcity. The 
difference between these two values of n is within the uncertainty quoted in sec. 3.2. 
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For the H resonance (Am^ ~ 3 ■ 10^^ eV^) and E = 10 MeV we have pres — 4000 g ■ cm~^, 
dres — 5 ■ 10^ Km (see fig. IT^ and — 5 s. The period Atg during which the H resonance is 
on the shock front (i.e. pi < pres < pf) can be estimated as the time which the shock wave 
takes to propagate from the position of the resonance, dres {pi = Pres) to the position df at 
which Pf = Pres- The latter can be found using the density profile ((Zj): 



Then the interval Ats equals: 



df = dresS}'^ . (118) 
At, = ^(^l/3_l)_ (119) 

Taking ,^ ^ 10 (eq. ()113j) ) and constant velocity Vs-, ()114j) . we get 

which is comparable with the arrival time ts- 

According to eqs. pi7|) and ()120|) . tg and At, increase with energy: for the lowest 
detectable energies in the spectrum one gets ts(6MeV) = 3 — 5 s, while for the highest 
energies ts(70MeV) = 7 — 10 s is obtained. 

Let us mark, however, that the details of the shock-wave propagation depend on the 
specific characteristics of the progenitor and therefore even large variations of parameters 
are expected between different SN models. 



Let us consider the conversion in the H resonance and its modifications due to the pas- 
sage of the shock wave. For a given neutrino energy three time intervals can be defined: 

1. The pre-shock phase, t{E) ^ ts{E) ~ 3 — 5 sec, in which the shock has not yet reached 
the resonance region and therefore has no effects on the neutrino conversion. 

2. The shock phase, t{E) = ts{E) — {ts{E) + Ats{E)) , when the resonance condition p = pres 
is fulfilled in the step of the shock front. 

As it follows from the fig. El in this phase three resonances appear at three different 
radii: two of them are on the walls of the hot bubble, while the third is at the shock front. 
We denote these resonances as Hi,H2, from the inner to the outer. 

For a given value of sin^6'i3, the conversion pattern due to the presence of these reso- 
nances exhibits a variety of possibile scenarios, depending on the adiabaticity character of 
these resonances, and therefore on the details of the density profile of the bubble and of the 
shock front. The shock-wave propagation has no effect if sin^ ^13^ 10~^, corresponding to 
non-adiabatic neutrino conversion in all the relevant resonances during both the pre-shock 
and the shock phases. Therefore in what follows we discuss the case in which sin^ ^13 is in 
the adiabatic region of the H resonance in the pre-shock regime, that is, sin^ ^13^ 10~^ (see 
fig. P), and the Hi and H2 resonances are adiabatic, while the adiabaticity is broken in the 
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resonance. In this case the permutation parameters p (for n.h.) and p (for i.h. ^) are 
given by eqs. (Pm) -()23 p with Ph being the jumping probabihty in the resonance. For 
maximal adiabaticity breaking 

~ ~ 1 • (121) 

The effect of the arrival of the shock-wave consists in a sudden change of the character of the 
H level crossing from adiabatic to maximally non-adiabatic, with significant modifications 
of the observed neutrino signal. 

3. The post-shock phase, t{E) > ts{E) + Ats{E), when the matter density at the shock 
front is smaller than the resonance density pres > Pf, so that only one level crossing, Hi, 
survives. The conversion effects are described by the transition probability in the Hi reso- 
nance: Ph — Phi ~ 0. In the transition between the shock and the post-shock phases, the 
resonance points H2 and H^ become closer until they merge and disappear ^. 

For the scenario discussed here we expect the following time dependence of the effective 
jumping probability Ph- 

'0 for t<tsiE) 
H ^ S 1 for t = ts{E) - (t, + At,{E)) (122) 
for t Z ts{E) + Ats{E) . 

In general, the shock wave effect is proportional to Pf/3, which in turn depends on the 
neutrino energy. Therefore not only the time but also the size of the effect depend on the 
energy. The ratio of the adiabaticity parameters of a resonance at and above the shock 
front (where the density profile equals that of the progenitor star) equals the ratio of the 
corresponding density gradients: 

7/_, _{dp/drh ..^o^ 
li {dp/dr)o 

From this and eqs. (|13I14|) it follows that for sin^ 6'i3 = 0.01 the adiabaticity is violated at 
all relevant energies if kgrad = 10^. That is the gradient of density in the shock wave is 4 
orders of magnitude larger than the gradient of the progenitor profile (at the same density). 
If kgrad = 10^, we get Ph^ ~ 1 for high energies: E ~ 70 MeV but Ph:^ ~ 0.5 for E ~ 6 

^The conversion in the non-resonant channel {De fon n.h. and for i.h.) can be modified by the 
shock-wave propagation if the density profile at the shock front is very steep so that the adiabaticity of the 
conversion is broken. However it can be checked that this requires unphysically large slopes of the shock 
front and therefore this possibility will not be discussed here. 

^Three low density (L) resonances may appear at these later times, i ~ 20 s, as shown in the figure 
1121 The one at the shock front could have a non-adiabatic character, so that the condition of adiabatic 
conversion (P^ = 0) adopted here is not valid in this case. However this happens in the latest part of 
the signal, when the neutrino luminosity is small. For this reason the effects of the Shockwave on the L 
resonance are not discussed in detail here. 
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MeV. For smaller gradient: kgrad = 10^, the adiabaticity is not violated for low energies, and 
therefore there is no shock wave effect, but it is still violated for high energies: Ph^ ~ 0.5 
for E 70 MeV. In general one expects larger effects for high energies. Results for other 
values of sin^ ^13 can be obtained immediately by rescaling eqs. ()1H|14|1 . 

7.2 Shock wave effects on the neutrino spectra 

Let us consider the influence of the shock wave on the neutrino spectra in more details. 
For a given energy E the effect of the shock wave consists in a change of the flux when 
the front of shock wave reaches the corresponding resonance layer with resonance density 
Pres{E) and then restoration of the original (undisturbed) flux after the front shock shifts to 
smaller densities. In what follows for simplicity we will assume that the density gradient in 
the shock wave is large enough so that the adiabaticity is strongly broken for all observable 
energies. 

In the case of normal mass hierarchy (H-resonance in the neutrino channel) during the 
interval Atg the shock wave leads to a change of the z/g flux exiting the star from completely 
to partially permuted: 

F°(E) ^ F,{E) = cos^ e,2F^{E) + sin^ OuF^iE). (124) 

Let us introduce critical energy E^ such that F^{E^ = F^{Ec). Then for E < Ec we have 
F^{E) < Fs{E), and for E > E^ the inequality F^{E) > Fs{E) holds. So, the pattern 
of perturbation of the whole spectrum can be described as follows: first the perturbation 
reaches the lowest energies of the spectrum and then propagates to high energies. We can 
say that the wave of perturbation of the spectrum propagates in the energy scale. It leads 
to an increase of the flux of the electron neutrinos according to eq. ()124j) when the wave 
propagates up to Ec- Then for E > E^. the pertubation decreases the flux according to ()124p . 
So, we can say that the shock-wave effects consist in a wave of softening of the spectrum 
which propagates from low to high energies (softening wave). When the wave reaches the 
highest energies of spectrum, the perturbation start to disappear at low energies, restoring 
the original spectrum: 

F,(E)^F0(E). (125) 

The spectrum becomes hard again and the restoration moves again from low to high energies. 
The antineutrino conversion will not be affected by the shock wave. 

In the case of inverted mass hierarchy the H-resonance is in the antineutrino channel 
and therefore the shock wave will modify the antineutrino conversion. Similarly to p24|) v^. 
spectrum changes as 

F^{E) ^ Fs{E) = sin^ 9i2F^{E) + cos' 9i2F^{E). (126) 

Although now the change of the permutation factor is stronger (from complete permutation 
to weak permutation: (1 — p) = 1 ^ sin' ^12), the effect on the spectrum can be similar 
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to the normal hierarchy case due to the smaller difference of F^{E) and F^{E) fluxes. The 
neutrino flux will not be affected by the shock wave. 

7.3 Shock wave and the Earth matter effect 

The shock wave propagation can influence the Earth matter effect. Indeed, for the normal 
mass hierarchy the Earth matter effect in z/g channel vanishes for adiabatic H transition 
while it is maximal for maximal violation of the adiabaticity in the H resonance, = 1 
13 inS] • The effect increases fast with energy, reaching 30 - 50 % at ~ 50 — 70 MeV. 

Taking this into account, for the scenario ()122p we get the following time dependence: 

in the pre-shock phase, no regeneration effect is realized in the z/g channel due to the 
adiabaticity character of the H resonance. The effect is small in the first few seconds of 
the shock wave phase: in this interval the adiabaticity is broken at low energies, where the 
Earth matter effect is small. After 5-7 sec the softening wave reaches the high energy 
region, where the Earth matter effect is large. So one expects the absence of the Earth 
matter effects in the first 5-7 sec of the neutrino burst and then its fast increase with time 
in the high energy part of the spectrum. We will call this the delayed Earth matter effect. 

Being unaffected by the H resonance, the Ue spectrum exhibits the Earth matter effects 
for the whole duration of the neutrino signal. 

In contrast, if the mass hierarchy is inverted, the Earth matter effect is observed in the 
neutrino channel during the whole burst and it appears in the antineutrino channel in the 
late stage (after 5-7 sec). 

7.4 Shock wave effects, mass hierarchy and ^13 

Precise measurements of the energy spectra in different moments of time can, in principle, 
reveal spectrum distortions which depend on time and propagate from low to high energies. 
For this however very high statistics is needed. 

Another possibility is to use some global characteristics of spectra and their change with 
time to look for shock wave effects. In particular, one expects 

• decrease of the average energy of the spectrum (E) which reflects its softening during 
the shock wave phase; 

• increase of the relative width F, as a consequence of the appearance of a composite 
spectrum during the time Atg] 

• appearance of the Earth matter effect in the late stage of the burst; 

• change in the total even rates [SB]- 

A problem of dealing with global characteristics is the possible degeneracy between 
the changes of the neutrino energy spectra due to shock wave effects and those due to 
astrophysical factors (cooling of the energy spectra, decay of the luminosity, etc.). 
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The shock-wave effects on the (or z/g) signal could be identified by studying the time 
dependence of the ratios of v^, (for n.h.) or z/g (for i.h.) CC to NC event rates at SNO 
Since the NC rate is not affected by shock-induced conversion effects, the CC to NC ratio 
represents a particularly "clean" quantity, in which various time-dependent features other 
than the shock effect (e.g. cooling of energy spectra) are at least partially subtracted. Using 
this method, in ref. jHSl it is found that the distortion due to the shock-wave amounts to 
~ 30%, whose statistical significance has not been clarified, however. 

An alternative approach to probe effects of the passage of the shock wave could be 
to study the time dependence of the ratio Rtot of the total number of Vf, and v^. events 
(eq. ()58j) ). where still an at least partial subtraction of cooling effects is realized. In this 
quantity the uncertainties related to the poor knowledge of astrophysical parameters can by 
no means mimic a specific time structure and therefore should be distinguishable from the 
shock effects. 

The observation of the delayed Earth matter effect is another signature of the shock 
wave effects. 

Let us summarize the possibilities to identify the mass hierarchy and measure (or restrict) 

^13- 

The observation of shock-effects in v^. (z/g) channel would select the normal (inverted) 
mass hierarchy and would tell that sin^ 6'i3 is relatively large, sin^ ^13^ 5-10~^, so that Ph < 1 
significantly in the pre-shock phase. Moreover, it would allow to study the physical features 
of the shock-propagation (speed, etc.). The exclusion of shock-effects would imply that 
sin^ ^13 is in the adiabaticity breaking region in the pre-shock times or that the properties of 
the shock wave differ significantly from the predictions (e.g. that the shock stalls, or travels 
with smaller velocity, so that the shock does not reach the resonance region during the 
duration of the burst, or that the shock front is not steep enough to change the adiabaticity 
character of the conversion) . 

Quantitative estimates of the bounds on sin^ Oi^ that can be obtained from the study of 
shock-induced time dependences of the signal will be given elsewhere |7Uj . 

8 Discussion and conclusions 

I. We have studied the effects of the 13-mixing on the conversion of the supernova neutrinos. 

If is negligible (sin^^i3^ 10~^), at the cooling stage a partial permutation of the 
original and u^, Ve and Vx energy spectra occurs in the star due to the large 12-mixing. 
For larger values of ^13, sin^^i3 > 10~^, these conversion effects are modified in the (z/g) 
channel for normal (inverted ) mass hierarchy. 

For normal mass hierarchy the degree of permutation in v^. channel increases from (1 — 
p) = cos^ Qyi ~ 0.7 — 0.8 to (1 — p) = cos^ 6^12 + (1 — Ph) sin^ 6*12, resulting in an hardening of 
the observed z/g spectrum. If sin^ ^13 > 10~^ the transition driven by dn is adiabatic and the 
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change is maximal: {1 — p) ~ 1. This corresponds to a complete permutation and therefore 
maximally hard u^, spectrum. 

Some additional distortion of the energy spectrum is expected for sin^6'i3 < 10""^ due to 
the dependence of the jumping probability Ph on energy. 

In the case of inverted mass hierarchy the 13 mixing leads to an increase of the de- 
gree of permutation in the antineutrino channel: from (1 — p) = sin^ 6*12 to {1 — p) = 
sin^ 6*12 + (1 — Ph) cos^ 6*12- Again for sin^ 6*13 > 10~^ the effect is maximal: (1 — p) ~ 1. 

In this paper we have elaborated methods to study the effects of the 613 mixing in pres- 
ence and absence of Earth matter effects. Various possibilities to probe the mass hierarchy 
and measure (or restrict) ^13 have been discussed as well. 

II. The permutation due to ^13 changes the average energy and width of the z/g (for 
n.h.) or z/g (for i.h.) spectrum. Signatures of these modifications can be obtained from the 
comparison of the features of neutrino and antineutrino spectra. 

In particular, for large 13-mixing and the normal mass hierarchy one expects {Eiy)) > 
{E{u)) and {T{u)) < {V{u)). Opposite inequalities are expected in the case of inverted 
hierarchy. 

In this connection we have studied the observed spectra of events from +d and Ve +P 
CC reactions at SNO and SK respectively; in particular, the following observables have been 
considered: the ratio of the average energies of the Ve and z/g events, te = (E) / {E), the ratio 
of the widths of the energy spectra, = T/F, and the ratio Rtaii of the numbers of events in 
the high energy tails (above certain energy cuts Ei and El) at the two experiments. These 
quantities depend on ^13 via the jumping probability in the H resonance, Puidn)- 

Although the 13-mixing produces a strong conversion, the corresponding observable ef- 
fects may not be so large due to presence of neutrino fluxes of all types and also due to 
effects of the 12-mixing effect (see item I ). Furthermore, substantial astrophysical uncer- 
tainties make the identification of the 13-mixing effects difficult, due to possible degeneracies 
between the parameters of the original neutrino spectra and the oscillation parameters. 

We have studied the effects of uncertainties on: (i) the energy spectra and luminosities 
of the neutrino fluxes originally produced inside the star, (ii) the density profile of the star, 
(iii) the mass square splitting Am|]^, (iii) the mixing angle 612- 

Methods to take these uncertainties into account in the data analysis have been elaborated. 

III. We have considered the cases of adiabatic transition {Ph = 0) with (A) of normal 
mass hierarchy and (B) inverted mass hierarchy and we have compared them with the case 
(C) of maximally non-adiabatic conversion (i.e. negligibly small Oi^: sin^ 6*13 < 10~^). 

The scatter plots of the observables {te, r^, Rtaii,---) for these three cases show that 
there are regions - in the space of these parameters - where only one or two possibilities 
exists, implying the exclusion of the remaining ones. 

In particular. 
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• The case A is excluded if observations give Rtau ^ 0.06 and/or rr ^ 1.1. If, in 
contrast, the experiments give RtaU > 0.22 the normal hierarchy would be identified. 
This result would be further supported if rr ^ 0.9 is also found. 

• The case B is excluded by (i) large values of RtaU, Rtaii ~ 0.07, (ii) large values of the 
ratio of average energies te ^ 1.2 and (iii) rr ^ 1. In contrast, the identification of 
this scenario appears difficult. 

• The case C is difficult to single out. The result rp ^ 1.1 and 0.06 ^ Rtau ^ 0.2 would 
exclude B and A, and therefore indicate that Ph > 0, corresponding to small values 
of ^13: sin^^ia^ few 10-\ 

IV. The high energy tail method allows, in principle, to establish the mass hierarchy 
and to restrict or to measure sin^6'i3. Taking the thresholds = 45 MeV and = 55 
MeV and tan^ 6*12= 0.38 as an example, we find that if i? > 0.22 the normal mass hierarchy 
can be identified and a lower bound on sin^ 6*13 can be obtained. For smaller values of the 
mixing the hierarchy can not be established; however strong bounds on sin^ ^13 are put if 
the mass hierarchy is know from other results (e.g. by measurements of other parameters 
of spectra, by studying the Earth regeneration effects etc.). A measurement of ^13 will be 
possible only if a reduction of the uncertainties on the original neutrino fiuxes is achieved. 

The Earth matter effect enhances the identification power of the method. Since the 
regeneration is larger in neutrino than in antineutrino channel, the sensitivity of RtaU to 
the n.h. (i.h.) is enhanced (reduced) with respect to the no Earth-crossing case. In partic- 
ular, in presence of Earth effects the method of analysis discussed in this paper would be 
complementary to the information that can be extracted from the observation of oscillatory 
distortions in the energy spectra of u^. and/or z/g events in the detectors. 

The method of the high energy tails can be improved if (i) additional information about 
parameters of the neutrino fluxes is obtained from the studies of the whole spectra (and not 
only of the high energy tails) (ii) if the dependence of the ratio Rtaii on the energy cuts is 
studied, (iii) higher energy cuts are used. These improvements can be implemented if the 
statistics is large enough. 

V. The study of shock- wave effects on the neutrino burst, which can be realized at 
late times (t ^ 5 s) could provide additional information on the 6*13 angle and the mass 
hierarchy (in support of the methods discussed in sections 5. and 6., which refers to early 
times, t ^ 5 s). In particular, the presence of time dependent spectral distortions produced 
by the shock-wave passage through the resonance layers would lead to the lower bound 
sin^ ^13^ few ■ 10~^ and establish the n.h. (i.h.) if observed in the Ue (t'e) channel. 

An important signature of the shock wave is the appearance of the delayed Earth matter 
effect in the neutrino or antineutrino channel at the late stages of the burst (t > 5 — 7 sec). 
The channel of appearance will identify the mass hierarchy. 
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Figure 1: The transition probability in the H resonance, Ph, as a function of sin^ 6*13 for two 
values of the neutrino energy (upper panel) and of the neutrino energy for different values 
of sin^^i3 (lower panel). We took C = 4 and jAmlgl = 3 ■ 10~^ eV^. The shaded regions in 
the upper panel represent the uncertainty associated to C = 1 — 15. 
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Figure 2: The survival probabilities p for n.h. (upper panel) and p for i.h. (lower panel) as 
functions of sin^ ^13 for tan^ 612 = 0.38 and different values of the neutrino energy. Other 
parameters as in fig. ^ The shaded region in the lower panel represents the error bar 
associated to a measurement p = 0. The corresponding lower bound on sin^6'i3, sin^^i3^ 
2 • 10~^, is small enough to be calculated in the zeroth order approximation (i.e. neglecting 
the explicit sin^ ^13 dependence in eq. fl23|l ). The error bar is much larger for the neutrino 
channel and normal hierarchy and was not drawn. 
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Figure 3: Scatter plot in the plane Rtaii — Rtot for the cases A, B and C discussed in sec. 
15.11 The quantity Rtot is the ratio of the total rates of ly^ and z/g events, eq. (jSHl)- The cuts 
El = 45 MeV and El = 55 MeV were taken for Rtau- 
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Figure 4: Same as fig. El in the Rtaii — S plane, where S is defined according to eq. 
We took the cuts E'j^ = E'j^ = 25 MeV for the calculation of S. 
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Figure 5: Scatter plots in the space of the variables: Rtaii, ^e, rr, eqs. (IK^ . (jUJ) and 
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Figure 6: The factor Q (region between the hnes), eq. (jHSI), as a function of the temperature 
Tx- The width of the region corresponds to the interval T^ — T^ = 0.35 -j- 0.5 MeV. The cuts 
El = 45 MeV and El = 55 MeV have been taken. 
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Figure 7: The dependence of the ratios a (upper panel) and a (lower panel), eq. on 
the energy cuts and respectively, for = 7 MeV, = = L^, rie = rj^ = rj^ = 
and different values of the z/g and Ue temperature (expressed in term of ratios in the legend). 
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Figure 8: The factor r, eq. (ff^ . as a function of Ph for tan^6'i2= 0.38 and both normal 
and inverted hierarchy (sohd hne). The bands within the dashed and dashed-dotted hues 
correspond to 10% and 20% uncertainty on tan^ 9i2 respectively. 
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Figure 9: The same as fig. |H1 where the assumption of hierarchy of the energy spectra has 
been relaxed to a = N^/N^^ < 1 and a = A^^/iV? < 1 (eqs. (IH21)-(IHSI)). The possible 
values of r are those within the solid contours; the dashed lines correspond to a = 0.17 
and a = 0.34. We took tan^6'i2= 0.38; for simplicity the effects of uncertainties on this 
parameter are not shown. 
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Figure 10: The ratio Rtaii as a function of sin^ 6*13 for normal mass hierarchy (sohd hues) and 
inverted hierarchy (dashed hues), in absence of Earth crossing. The bands within the thick 
hues represent the values that Rtaii can take once all the possible theoretical uncertainties 
are considered (see text). The thin lines refer to the "ideal" case a = a = 0. We have taken 
tan^^i2= 0.38, = 45 MeV and = 55 MeV. Three examples of experimental results 
for Rtaii are shown with two different (98% CL.) error bars (marked by bullets and stars) 
corresponding to T^, = 7 MeV, L^, = Lg = = = 5 ■ 10^^ ergs, and D = 8.5 kpc and 
D = 4 kpc respectively. 
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Figure 11: The same as fig. Elfor Earth crossing trajectories. We took Amgi = 5-10 ^ eV^ 
and the nadir angles On = 84.3° at SNO and 9n = 24.5° at SK. 
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Figure 12: The density profile at various times post-bounce as adopted in ref. [22] • The 
horizontal lines represent the resonance densities pH (for two values of the neutrino energy) 
and pl'-i the three high-density resonances Hi,H2,H3 are marked for t = 5 s and E = 10 
MeV. 
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